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Abstract 

We present and analyze a penalization method wich extends the the method of [2] to 
the case of a rigid body moving freely in an incompressible fluid. The fluid-solid system is 
viewed as a single variable density flow with an interface captured by a level set method. The 
solid velocity is computed by averaging at avery time the flow velocity in the solid phase. 
This velocity is used to penalize the flow velocity at the fluid-solid interface and to move 
the interface. Numerical illustrations are provided to illustrate our convergence result. A 
discussion of our result in the light of existing existence results is also given. 

>" 

00 

^ ■ 1 Introduction 

In this paper we are concerned with the numerical analysis of a penalization method for the two- 
way interaction of a rigid body with an incompressible fluid in three dimensions. The traditional 
numerical approach to deal with fluid-structure problems is the so-called ALE (for Arbitrary 
Lagrangian Eulerian) method where fluids (resp solids) are described in an Eulcrian (resp La- 
grangian) framework. Fluids are computed on a moving mesh fitting the solids and stress and 
velocity continuity are used to derive the appropriate boundary conditions on the fluid/solid in- 
terface. A convergence proof of a finite-element method based on this approach can be found in 

Alternate methods can be devised where the whole fluid-solid system is seen as a multiphase 
flow and the fluid/solid interface is captured implicitly rather than explicitly. Likewise, the inter- 
face continuity conditions are recovered in an implicit fashion. The rigid motion inside the solid 
phase can be enforced through a Lagrange multiplier [TT] . The method we consider here is of this 
type but the rigid motion is approximately satisfied in the solid through penalization. 

Penalization methods have already been considered in the past for this problem. In [6] the 
authors considered a single solid ball and worked inside the frame moving with its center. Then 
they penalized the mean velocity of the (virtual) fluid inside this ball. Their method is restricted 
to one ball. In[Tni[T5] the penalization is applied to the deformation tensor inside the body. In [THj 
this method is used to prove the existence of solutions for the fluid-solid interaction variational 
problem in two dimensions. In [15j it is used together with a two-dimensional finite element 
method in a variational framework. Here the penalization is applied to the flow velocity itself. 
The method thus extends the one devised and analyzed in [5] in the case of a rigid solid with 
prescribed motion. 

In our method the determination of the body velocity is part of the problem. This velocity, 
instead of the flow velocity, is used to move the solid phase. This has a crucial practical importance, 
in particular for problems with large displacements and strong shear, since it ensures that the 
solid remains rigid at the discrete level, although the rigidity constraint in the flow field is only 
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approximately satisfied. A vorticity formulation of the method and its validation on a number of 
2D and 3D reference cases are given in [7J. An outline of the paper is as follows. In section 2 we 
recall the weak formulation of the problem and we describe the penalization method. Section 3 is 
devoted to the convergence proof. In section 4 we provide some numerical illustrations. Section 5 
is devoted to some concluding remarks. The proofs of some technical results used in section 3 are 
given in the appendix. 

2 Weak formulation and penalized problem 

Let fl be an open bounded domain of K 3 , filled with a viscous incompressible and homogeneous 
fluid of density pf > and viscosity p, > 0. Inside this domain, we consider the motion of an 
immersed homogeneous rigid solid of density p s > during a time interval [0, T], T > 0, chosen 
so that the solid never comes in contact with dil. For £ S [0, T], we denote by f2/(t) and Cl s (t) the 
non-empty fluid and solid open connected domains, with f2 s (£)Ufi/(t) = f2 and f2y(£)nfi s (£) = 0. 
The center of mass of the solid is denoted by xg(£), its mass and inertia tensor by M and </(£). 
Without loss of generality we assume that M = 1. Then 



£g(£) = / p s xdx, J(t) = / p s (r I — r <g> r) dx 
where r(x, £) = x — xc(t). The system is subject to a body density force g (usually gravity). 
2.1 Weak formulation 

The basic formulation of this fluid-solid coupling is the following : given initial conditions, 

x' G (0)=v° g , w u (0)= W °, u = u°, fi 5 (0)=fi° (1) 

supplemented with 

x G (0)= xdx, X s (x,0)=x, (2) 

find £ — > n s (t) and (x,t) — > (u(x,£),p(x,£)) solution for £ > of 

p/(w t + (u • V)tt) - 2/xdiv(D(u)) + Vp = on %(£), (3) 

divu = onfi/(i), (4) 

u = on SO, (5) 

ti = ig + w„ x r on 9fi s (£), (6) 

*g(*)=S+ / (Sn)ds, (7) 
Jdn B {t) 

J(t)u' u (t) =-w„(t) X (J(t)w tt (t))+ / p s (rxg)dx+ I rx(Sn)ds, (8) 

Jn a (t) Jdn s (t) 

n s (t) = x s (t,n°), (9) 

9A 

■^= I ' G (t)+4)xr(I I (£),£), (10) 

where n denotes the unit outward normal on <9f2 s (£), and S is the fluid stress tensor. 
In this formulation the last two equations describe the rigid motion of Q s (£). In order to give a 
weak formulation of this problem, let us introduce some function spaces. From now on, u will 
denote the velocity field on the whole computational domain f2. We define 

V = {u e H%(n), divit = 0}, H = {u e L 2 (Q), divu = 0, u< n = on dtt} 
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and, with the notations of [TH] extended to the three dimensional case, 

K{t) ={ii£V, D{u) = in Q s (t)} 

= {ueV, 3{Vu,uj u ) £ I 3 x 1, u = V u + uj u x r in fl a (t)}. 

Next we define a density on the whole domain by setting p — p s Xn s (t) + PfXn s (t), where \A 
denotes the characteristic function of set A, which takes value 1 inside A and outside. Let 
us note Q = f2x]0, T[. Then the weak formulation is the following [14]: given initial conditions 
H° = Xn°, P° = P s H°+pf(l-H°) and ii = ii°e /C(0), find (x,t) -► (p(x, t), u(x, t), H(x, t)) such 
that 

'«£ i°°(0,T,H)n J L 2 (0,T,V), H,peC(Q,T;L"(n))Vq> 1, 
u(t) G K(t) for a.e. t e]0,T[, with Sl 3 (t) = {x G O, = 1}, 

V£ e fl-^Q) n £ 2 (0,T;£(t)), 



[pit • <9 t £ + (p(u ■ V)u - 2pD(u)) : £>(£) + pg ■ £] dx 



d 
dt 



pit • £ dx, 



WgC^Q), ^(T) = o, 

H^- + Hu- Vibdxdt 
at 

dib 

p— — V pu ■ Vip dxdt + 
o Jn 9t 



(11) 



H°ip{0)dx = 0, 
pV(Q) da; = 0. 



Note that we could equivalently have defined p = p s H + pf(l — H), as p° is piecewise constant, 
and transported by the same velocity field than H. 



2.2 Penalized problem 

For r\ > 0, we consider the following penalized problem: 

given (p n (0) = p°,u„(0) = u°,H n (0) = Xu°)> to find (p v , u v , p v , H v ) , with 

p^Hr, Gi°°(]0,T[xO), u IJ G£ oo (0,r;W)nL 2 (0 ) T;V), p„ e L 2 (Q) 
solution on Q of 

/ Q u \ 1 

Pv \ + ( u n ' V K; J ~ 2fj,dxv(D(u v )) + Vp v + -p n H n {u n - u, hS ) = p v g 



divw r) = 



Prjir^ x u n )H n dx x r n 



Pvt 



vt 



+ u-n-Vpn = 

U^g.VHr, = 



(12) 

(13) 

(14) 

(15) 
(16) 



We set fi^(i) = {x G fi, H v (x,t) — 1}. In equation (jT5J) we divided the first term by M v = 
Jn p v H v dx, which is not constant in time in general. On contrary we have |f2y(f)| = J n H v dx = 
since u v )S is divergence free and if,, vanishes on 917 (we assumed no contact of the solid with 
dft). The inertia tensor is defined as 



with r r , 



x - x Gl) 



p^H^r 2 ! - r v <g> r v ) dx 
In Pv H v x dx - 



p,,(r 2 I — r,, ® r,,) dx. 



For a G M 3 \{0}, a T J„a — } n vn\ pr^fn x a\ 2 dx > min(p s ,p/) J^>?/^ x a| 2 dx (see estimate ([T8])). 
This last quantity being strictly positive for an open nonempty integration set, is nonsingular 
(we recall that = > 0). 
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Before stating our convergence result, a few remarks are in order. 

First one may wonder about the well-posedness of the above problem. However it will directly 
result from the a priori estimates and convergence arguments given in the following that this 
problem does have at least a weak solution. Indeed, these arguments could easily be used to show 
the convergence of the solutions to a linearized version - or finite-dimension approximation - of 
(fT"2 | -(|16 p . Next we can observe that in this model we penalize the difference between u v and the 
projection of onto velocity fields rigid in the solid domain, namely u n ^ s (see lemma l3Tl below) . 
The density is transported with the original velocity field so that estimates on the Navier-Stokes 
equations are easier to obtain. The characteristic function is transported by the rigid velocity 
so that the shape of f2y(i) remains undeformed (this is exactly the Eulerian counterpart of (J5J- 
UOp V As observed in [7j this has a practical importance (in particular it means that the rigid 
solid can be recovered exactly through simple algebra from its initial shape). As far as numerical 
analysis is concerned, it also provides " for free" regularity properties on the computed rigid body, 
as soon as the initial body is smooth. The price to pay is that the level sets of p v and H v do not 
coincide, i.e. in general we do not have p n H n — p s as in the non penalized formulation. Note also 
that in principle we should prescribe a boundary value for on dQ when u n ^ s is inward. Since 
our analysis is restricted to times when the solid body does not approach the boundary of the 
computational box, we can take this boundary value to be zero, which amounts to solve (|16| on 
W 1 and take its restriction to fi. 

In the following sections we will prove the convergence of at least a subsequence of {p v , , p n , H n ) 
to the weak solution defined above. Next section starts with some a priori estimates which will 
provide weak convergence of subsequences. In section 1331 we will have to use more sophisticated 
tools adapted from [TS] to get some strong convergence in u v which will allow us to pass to the 
limit in nonlinear terms of {P^)- More precisely we prove the following result. 

Theorem 2.1. Under the regularity assumptions of section^ let (p v , u^^p^, H^) a solution of 
(Pr/)- Then there exists a subsequence of (p^,^, H^) and functions (p,u,H) such that 

Prf — ► p, H n — > H strongly in C(0, T; L q (£l)) for all q>l, 

u v -> u strongly in L 2 (Q) and weakly in L 2 (0, T; Hq(Q,)) n L°°(0, T; L 2 (il)) 
and such that (p,u,H), is a solution of 

Before proceeding to the proof, let us point out a few remarks. For a sake of simplicity in the 
notations we have stated our penalization method and theorem for a single rigid body. It will be 
apparent from the proof below that it readily extends to the case of several bodies. Furthermore, 
the time to which the convergence result is restricted, is essentially the time for which contact of 
the rigid body do not touch the boundary of il (in the case of several bodies it would be the time 
on which we can ensure that contact between bodies do not happen). As a result if we consider 
periodic boundary conditions and a single body convergence holds for all times. 



3 Proof of theorem 12.11 

The following lemma states that u^ iS , as defined in (P v ), is the projection of u v onto velocity fields 
which are rigid on Q^(t). 

Lemma 3.1. Let £ be a rigid velocity field, i.e. such that £(x) = +W{ x r(x) for some constant 
vectors € M 3 and £ M 3 . Then if u n , s is defined by \14\) there holds 

I PrjH^Ujj - u v . s ) ■ £_dx = 0. (17) 

Jn 

Moreover, the result holds if £ is a time dependent velocity field rigid in f2J?(t) at time t. 
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Proof. Let the mean translation and angular velocities be defined as 

V u = — I p v H v u v dx ui u = J^ 1 I p v H n (r n x u n ) dx 
M n Jn Jn 

then 

p^H^u^ - u rhS ) • £cfe = I p v H ri [u v - (V u + uj u x r v )] ■ [Vj + U£ x r v ] dx 

Jn 

= Vj ■ / p v H n u v dx + lo^ ■ / p v H v (r v X. u v ) dx - V u - Vjt / p v H v dx 
Jn Jn Jn 

- Vu- w s x / p v H v r n dx) -V ( ■ [uj u x / p^H^r^ t 

\ Jn J \ Jn 

- / p v H n {u u x r n ) • x rv,) da; 
Jn 

= V e ■ (M V V U ) + w £ • (J, w„) - K u • (M„ V 5 ) 

- V u • ^ x J prjH^ dx^j - Vt ■ (u u x J p v H v r n dx 

- / p n H n {u u x r. n ) ■ (u£ x r n ) dx. 
Jn 

As (w„ x r n ) ■ (uj£ x r,,) = (w^ • w„)r^ - (r v ■ uJ^)(r v ■ uj u ), we have / p v H v (uj u x r v ) ■ (w 5 x ry,) dx = 
Finally, by definition of r v , / p n H v r n dx = 0, and we get 

p v H v (u v - u v<s ) ■ £dx = W£ • ( J n lj u ) - W£ ■ ( J v oj u ) = 0. 

n 

□ 

3.1 Estimates for transport and Navier-Stokes equations 

In all the sequel, C denotes a positive constant. At this stage, we consider a given time interval 
[0, T]. The value to which T must be restricted will be given later in this section. 

Standard estimates for transport equations (|15j) and (|16j) show that p v and H n are bounded 
in L°°(0, T; L°°(fi)). More precisely, for all time t G [0, T] , 

p m in ■= mm(p s , pf) < p v (x, t) < max( j o s , pf) H n (x, t) € {0, 1} a.e. ie!l. (18) 

Thus, up to extracting a subsequence, we can assume that 

p v -± p in i°°(0, T, L°°(Cl)) weak*, (19) 

and 

H n ^ H in L°°(0, T, L°°{Cl)) weak*, (20) 

where and p satisfy the bounds (fT8|) . We set H s (t) = {i £ fi, H(x,t) = 1}. Concerning 
Navier-Stokes equations, multiplying (| 12(1 by it,, and integrating on Ct, we get: 

/ Pv ( "T"^ ' ) ' u n dx — 2p ( dxv(D(uri)) ■ u v dx + f it,, ■ Vp,, dx 

Jn \ ot ' J Jn Jn 

+ — / PriH v (u v — u v , s ) ■ u v dx — I p v g ■ u v dx. 
V Jn Jn 
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Classically we have from incompressibility and homogeneous boundary conditions on u v , 

L Pv + K ■ VK ) ■ Ur > dx=l *L d[Pv dt ?) dx - (2i) 

From Lemma 13. II we get 

: Jn 
and since u v is divergence free and vanishes on dfl, 

/ u v ■ \7p v dx = 0. 
Jn 

Collecting terms we get, since from (fT8|) y/H v = H^, 

1 d 2 2 1 2 - 

2^llv^^lll2 ( ^ ) +MPK)ll L 2 (0) + -|lV^' J ff I) (M I) -^, S )|| L2(0) < || v / ^"» 7 IU 2 (n)ll5llL°°(Q)llp»jlli2 ( Q ) 
which upon time integration on [0, T] gives 

J ;,o, "■■ -/'Hi-'!",,) , "i- 

V 

2 



2 

+ 2fx\\D(ur,)\\ L2(Q) + -\\^H v (u v -u V!a )\\ L2{Q) 



, r T — 

< \\VPvo u vo\\L2(n) + C \\VPn( s ) u v( s )h 2 (n)ds. 



Applying Gronwall Lemma, Poincare inequality and bounds from (|18[) gives the following esti- 
mates : 

u v bounded in L 2 (0, T, Hq(Q)), (22) 
V/V^ and Uyj bounded in L°°(0, T, L 2 (fi)), (23) 

—y/p H v (u n - U v , s ) and -^=H v (u v - u^ a ) bounded in L 2 (0, T, L 2 (f2)). (24) 
Thus we can extract subsequences from p n , and H^, still denoted by p v , u n and i/^, such that 

— ^ u in i 2 (0, T, flo( n )) weak ' ( 25 ) 
Tp^u,, ->> x and u, n u in i°°(0, T, L 2 (fi)) weak*, (26) 
yfp^H v u v - yfp^H^Ur,^ -* and fl^it,, - H v u v , s -> in L 2 (0, T, £ 2 (f2)) strong. (27) 

The identification of x with y/pu results from strong convergence results proved by Lions and 
DiPerna on transport equations. [9] theorem II. 4, (|25p and incompressibility imply 

p n -> p in C(0, T, L 9 (f>)) strong Vg G [1, +oo[ (28) 

with p solution of 

jp t +u-S7p = Q onOx]0,T[, 
\p = po on £1 x {0}. 

From this strong convergence we can pass to the limit in the product ^fp^u^ : given v £ 
L q (0,T;L r (Cl)) with q > 2 and r > f, we write 



(■y/p^Uri — yfpu)v dxdt = / (u r; — u)y/pv dxdt + / (^/p,, — y/p)u v v dxdt. 
o Jn Jo Jn Jo Jn 

From the injection of H 1 into L 6 in dimension less or equal to 3 the first integral converges toward 
0. For the second integral we use the strong convergence (following ([25])) of Jp~^ in L s for a s such 
that u v v is in L s where s' is the conjugate exponent of s. Thus we have 

p^Ur, in i 9 (0, T, L r (n)) weak, for all g < 2, r < 6. (29) 
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3.2 Setting T and passing to the limit in the rigid velocity 

This rigid velocity is defined by 

Ur,,s( x > *) = u n,G (*) + u v (t) x r v {x, t), 

with 

u n .a(t) = — / p v u v H n dx and uj ri (t) = J" 1 / p v (r ri x u^H^dx. 

First we note that M v — J Q p n H v dx is bounded from below independently of 77 since p n is bounded 
from below and J Q H v dx = > does not depend on 77 or t. From the bounds on p v , H v and 
u n it is straightforward to show that 

u v . G {t) bounded in L°°(0,T). 

Likewise, from the definition of J v we observe that for a G R 3 \ {0} 



a T J n a > min(p s , p/) / \r, q x a\ 2 dx > 0. 



Moreover, the initial solid is regular and transported by a rigid velocity. We thus know that there 
is a ball of radius R > centered on the center of gravity xq^ included into ^(t). Then the above 
estimates implies 

a T J n a > min(p s , pf) / \r v x a\ 2 dx — min(p s , pf) / \x x a\ 2 dx = C(R)\a\ 2 

JB(xa„,R) Jb(0,R) 



with C{R) = > 0. Taking a = J n 2 b{J v is symmetric) we get for all tel 3 \ {0}, 

b T j-^b = \j-h\ 2 <^ r) \b\ 2 , 

which proves that each coefficient of J^ 1 is bounded independently of 77 and t. From the bounds 
on u v , H v and p n this implies that 

LJ v (t) is bounded in L°°(0,T). 

In particular this implies that the solid velocity u VtS is bounded in L°° by some constant M 
independent of r] and time. We can now define the maximum time for which the convergence 
result will be proved. If we denote by do the initial distance between solid and the boundary dSl 
then choosing for instance T = d /2M ensures that the body will not touch the boundary for 
t € [0, T]. In all the sequel we will assume this value of T. 

From the above estimates we can ensure that there exists uc(t) and u)(t) in L°°(0, T) such 
that, up to the extraction of subsequences, 

u v , s u s := u G +uj x r in L°°(0, T, weak*. 

Now we point out that taking the gradient of the rigid velocity field u VtS gives 

-ujI 

Vu ViS 

so that the Vm^ iS (and all subsequent space derivatives) is also bounded in L°°(0, T; L°°(Q)). In 
particular 

u v . s -± u s in L 2 (0, T, W hoc (Q)) weak*. (30) 
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We now wish to prove that u s (or equivalently uq and ui) has a similar structure as u^.s, 
that is, to pass to the limit in the expression defining u v ^ s . Using (|30p . the already mentioned 
compactness results of [5], applied to the transport equation on H v now gives 

H v ^> H a.e. and in C(0, T, L P (Q)) strong Vp € [l,+oo[ (31) 

with verifying 

jH t + u s .VH = on E" x (0,T) 
1 i? = H on E™ x {0}. 

Note that this Cauchy problem has been set in E™ because u s does not vanish on dtt, but H 
vanishes outside St. However we can prove by passing to the limit in ([27]) that Hu = Hu s . Thus 
div(tt-ff) = div(u s H) and -ff verifies a transport equation with velocity field u on f2 (note that no 
boundary conditions are needed on <90 since u is zero on the boundary). 

This convergence gives us the strong convergence of r v in C(0, T, L p (fl)), \/p > 1, and from 
|). (|31|) and (|25j) we can easily pass to the limit in the expression of u„ g an d w 7) to get 

puH dx 
pH dx 



; pun ax , . -i 

■"(?(*) = ^7 and w(t) = ( p{r 2 l - r (g> r)H dx) / p(rxu)Hdx. 



3.3 Strong convergence of 

The remaining part of the proof is more technical since we aim to prove the strong convergence 
of at least a subsequence of u v in order to be able to pass to the limit in the inertial term of 
Navier-Stokes equations. Classically, this is obtained thanks to a Fourier transform in time which 
provides an estimate on some fractional time derivative of u v which brings compactness [161 122] . 
Here these technics can not be used since the solid is moving. We instead rely on tools developed 
in [19]. 

Thereafter we will use, for a > and r S [0, 1], the following notations 

• n Sia (t) = {xe n,d(x,n s (t)) < a], 

• V° = {v e L 2 (Q), divv = 0, v ■ n = on dtt}, 

• V = {ve H r (ti), divv = 0, u = 0on dfl} (for r > 0), 

• K, r a (t) = {v(t) e V, D{v(t)) = in V'{Q s ^(t))} (which is a closed subset of V r ), 

• P^{t) the orthogonal projection of V on K. r a (t). 

To prove the strong convergence of a subsequence of u v in L 2 (Q) we write 

v — u\ 2 dxdt < [ / ( \p(u 2 — u 2 )\ dxdt + j f \2pu ■ (u — u v )\ dxdt ) . 

Jfl Pmin yJo J SI JO JO J 

From ([26]) the second integral on the right side converges to 0, thus 
r T r i / r T r <"T 

|2 7 j. ^ 1 ( / / I „ „,2 „„,2 



/ / \ u n ~~ u \ dxdt < / / \p, 1 u rj — pu | dxdt + / \(p n — p)u v \ dxdt . 

>0 Jn Pmin \Jq Jn J0 Jfl ' I 



S 



Moreover, by (I22p and (f2"5| the second integral on the right hand side converges to 0, and 



T 

Jn ' Pmin \Jq Jn 

T 



u v — u\ 2 dxdt < I / / \p v u v ■ P^(w^) — pu ■ P^{u)\ dxdt 



< 



\prjUr, ■ (u v — P^.(u 71 ))\ dx dt + / / \pu ■ (P£ (u) — u) I dx dt + l v 
a Jn Jo Jn . 

7-dlW ' K( u v) - P u ■ K{u)\\ L i{Q) 
Pmin 

+ \\pr,\\L~(Q)\\u v \\Li(Q)\\Pa( u n) ~ U v\\^(Q) 

+ \\p\\l<*>(q)\\u\\l2(q)\\PZ(u) - u\\ L 2 {Q) +l n 
where l n — > when i] — > 0. Finally, as p n is bounded in L°°(0, T, L°°(£i)) and using (f2"2"| . we get 

T f 1 

/ |u„ - u\ 2 dxdt < (\\p v u v ■ - pu ■ pr(u)\\ L i {Q) 

Jn Pmin 

+ C\\ P ai u r,) ~ V-Jl^Q) 

+ C\\P:{u)-u\\ L 2 [Q) +l ri (32) 

This decomposition shows that the sought convergence essentially amounts to prove that (up to 
the extraction of subsequences) 



]im\\PZ(u) - u\\ L 2 (Q) =0 (33) 



<T^0 

lim lim - NlU 2 (Q) = (34) 

lim lim • P^w,,) - pu ■ P^(u)\\ L i {Q) = (35) 

a— *0 77— >0 

To prove (|33H35[) we will make use of some Lemma that we state now. 

Lemma 3.2. Let (/„) be a sequence of functions bounded in L p (0,T) for some p > 2 and con- 
verging to almost everywhere on [0,T]. Then (f n ) converges strongly to in L 2 (0,T). 

Proof. Let e > 0. From Egorov theorem ([5], p. 75), the almost everywhere convergence implies 
that there exists A £ c [0, T] such that 

(\[0,T}\A e \ < e 

I /n ~~ * uniformly on A E 



which means 
Therefore 



3N G N,Vn > JV.Vt G |/„(£)| 2 < e. 
/ (/n(*)) 2 * < el^l < eT. 



As from the bound in L p norm, 



(f n (t)) 2 dt < / l"dt / {f n {t)fdt < e^C 

[0,T]\Ae J[0,T]\Ae J[0,T] 



with i + i = |, we get 



T 

2, 



{f n {t)Ydt<eT + e q C 

which proves the L 2 convergence. □ 
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Lemma 3.3. Letu(t) S H 1 (il s , cr (tj) such that W|gn» „(*)(*) = g(t) and (w(t),p(t)) € H 1 (il\il s , cr (t))x 
L 2 (fl\fl Sir7 (t)) solution of the Stokes problem 



-Aw(t) + Vp(t) = on n\n s<(T (t), 

div w(t) = on Cl\fl s ^(t), 

w(t) — git) on 8Q Si(T (t), 

^w(t) = on d$l. 

Then there exists Co > and C > 0, such that for all a < <jq, we have the following estimate: 

IK*)lli=(n\n^(t)) < c 'll«WllL2 ( o s , <T (t))l|VM(t)|| L 2 ( ^ At)) 
The proof of this result is postponed to the appendix. 
Lemma 3.4. The limits H , u and u s defined in \31\), V25\) and iS0\) verify 

Hu = Hu s . (36) 

Proof. First we claim that 

H v u v -± Hu in L q (Q, T, L r (il)) weak, with q < 2 and r < 6 (37) 
Indeed, let us introduce v £ L q (0, T, L r (Q,)) with q > 2 and r > |. We have 

/ / (H v u n — Hu) ■ v dxdt = I / H(u v — u) ■ v dxdt + I / (H v — H)u v ■ v dxdt 
Jo Jn Jo Jn Jo Jn 

From the injection of H 1 into L 6 in dimension less or equal to 3 and (|23|) we get 

u n u in L 2 (0,T,L 6 (n)) weak, 
and, since H is bounded in L°°(Q), 

lim / / H(u n — u] ■ v dx dt = 0. 

Moreover from (|3Tj) we easily get 

lim / / — H)u„ ■ v dx dt = 0. 

We next show that 

H v u VtS Hu s in V(0, T, L p (n)) weak, Vp € [1, +oo[. (38) 
Let v € LP(0, T, L p (Cl)) with p > 1. We write 

/ / (H^Urf^s — Hu s ) ■ v dxdt = I / H{u v , s — u s ) ■ v dxdt + / (H n — H)u v , s ■ vdxt 
Jo J si Jo Jn Jo J si 

As H is bounded in L°°(Q), with ((30]) we get: 

lim / / H(u v s — u s ) ■ v dx dt = 
v-^Jo J SI 

In addition, from (f3T|) we get 



lim / / — H)u„ s ■ v dx dt = 0. 
'J- Jo Jn 
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By (JSTD and (38]) we thus deduce 

H v u v - i^u,^ Hu- Hu s in L 9 (0, T, L r (fl)) weak, with g < 2 and r < 6 (39) 
Finally we recall that 

H v u v — H^u^ g — ► in L 2 {Q) strong, andthedesiredresultfollows (40) 
The result is finally obtained by identifying the limits in ((39)) et (|40[) . 

□ 

Lemma 3.5. 

V(t>0, 3^0 >0, V?7<r/o, Vte[0,T], fij(t) C Si(7 (t) and Q s (t) C nj j<r (t). (41) 
Proof. From pip with p = lwe have 

Ve > 0, 3770 > 0, Vr/ < Vte[0,T], \H v (x,t) - H(x,t)\dx < e, (42) 

Jo 

which means 

Ve>0, 3r/ >0, V77 < tjq, Vte[0,T], |fi*(t) \ fi.(t)| + \«Z(t)| < e. (43) 

By contradiction we suppose that we can find o~o > such as V770 > 0, there exists r\ < ijq and 
( e [0, T] for which at least one of the inclusions of (l4Tj) is false. This means that we can find 
x v (t) 6 fig(t) sucn as d(x v (t),fl s (t)) > a . fl^(t) is a rigid deformation of £1!?(0) so its boundary 
is C 2 . Thus, there exists a sufficiently small p independent of 77, such as for each point of fij(t) 
there exists a ball of radius p > containing this point and included in Sl!?(t). Then there exists 
also a ball of radius p :— min(p, <to/3) containing the point and included in O^(t). This latter ball 
is included in fi^(t) \ fl s (t). Indeed it contains a point at distance more than Go from £l s (t) and 
its diameter is less than 2cto/3. We thus got that 

3a , V770 > 0, 3t7 > 0, 3t e [0,T], |flj(t) \ n,(t)| > np 2 

with p independent of r] and t. This contradicts (|4"3"|) . □ 

Lemma 3.6. 

lim ||P;(u) - u|| £3( o,r,vO - 0, Vr e [0, 1[ (44) 
The proof of this result is postponed to the appendix. 

The next lemma is essentially, which is also proved in the appendix, is a rephrasing of the 
previous one with u v instead of u. The difference is that we do not have anymore u v — u^.s = 
ins ft%(t), but we do have an estimate on it, from which allows to pass to the limit as 77 goes 
to 0. 

Lemma 3.7. 

lim lim ||P;K) - zgU^o.T.v) = 0, Vr € [0, 1[ (45) 

a— >0 r\— >0 

Lemma 3.8. 

lim lim \\prjUr, ■ P^(u v )- pu- P^(u)\\ L im T ]xa) =0, We]0,l[. (46) 
Proof. Let r g]0, 1[ and er > 0. From Lemma [3~5l there exists 770 > such that V77 < 770, 

^(t)cf! v/3 (t), Vte [0,T]. 
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Let T) < t]q. Arguing as in [T^ we split [0, T] in Nt subintervals Ik = [(k — l)r, kr], r = T/Nt, 
k = 1, .., Nt- We choose Nt large enough (depending on a) such that 



n.,«r/2(frr) C n. )0 .(t) and n, i<r/3 (t) C ^, ff / 2 (fcr), Vt G J fc , Vfc = 1, ..,N T . 



(47) 



This is possible since Cl Sta (t) is moving with a rigid velocity field, L 2 in time. Thus the flow 
generated by this velocity field is continuous in time, and f2 SjCr (t) is the image of a through this 
map. On each subinterval Ik, k = 1, .., Nt, we consider the momentum equation 

du 1 
/°n"7^ + Pri(uri ■ V)u„ - 2^div(D(u J? )) + Vp v + -p v H v (u n - u rhS ) - p n g = 0. 

Let us consider a test function £ vanishing outside Ik and such that £(., i) € ^i/2(^ T ) f° r ^ ^fe- 
Since £(.,t) is rigid on n S)0 ./ 2 (A;T) D O s ,<t/3(*) ^ ^2(*)> Lemma [XT] yields 



[Pt?^ ' St + ® "77 - 2(iD(u v )) : D(f) + p^g • f] dxdt = 0. 

and (f2"5| we derive the following estimates: 



'h Jn 

From bounds given by |T 



/ D(u n ) : D(£)dxdt 
Ji k Jo, 



< \\D(u ri )\\ L 2 {IkL 2 {n)) \\D(0\\mi k ,L^n)) 

< C U\\L*(I k ,lCl /2 (kT)), 



and 



{prjUrj <g> u n ) : D(£) dxdt 



/ / Pv9 ' £ dxdt 
Ji k Jn 

Collecting terms we get 



< 



< 



< 



IKIU°°(fi)IK ® u v\\L 2 (n)\\ D (Q\\L2(n)dt 



\\pv\\L°°(I kl L°°(fl)) / \\Ur,\\ L 4(Q)U\\Hl(a)dt 



c / \\ u v\\h { n)\\ Vu rt\\h { Q)\\^\\H^n)dt 



ik 



< CIKIIJUd^cQ)) / l|Vu,||£ a(n) ||Cllir i(n)dt 



< C\\Vu ri \\l 2{IkL2(n)) \\i\\ L 4 {Ik , H l {a)) 



< 



C K\\L*(I k ,Htj(Q)) < C U\\mi k Xl /2 (kr))^ 



< \\Pv\\L~{I k X<~(Q.))\\i\\mi k Xl /2 {kT)) < C||fllL4(/ fc ,/Ci /2 (fcr))- 



i k Jn 



p v u v ■ £,t dxdt 



< C M\\Li(I k ,Kl /2 (kT))- 



As £(.,*) € K} aj2 (kr), &(-,*) € ^ /2 (fcr) C /C° /2 (fcr) and we have 



<r/2V 



IW^>&}I = 1^^,^/2(^)6)1 
= \(P"/2(kr)( Pv u v ),^}\ 

= \{j t P° /2 (kT)(Pr,U v ),0\. 
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Therefore 



which means that 

j t P° /2 (kT)( Pr , Uri ) bounded in L*(/ fc) (JC^kr))*). (48) 
Moreover p v u v is bounded in L 2 {Ik 1 L 2 (Sl)), 

P°/2(kr){p v Ur,) is bounded in L 2 (I k , /C° /2 (fcr)). (49) 

Since /C° / 2 (fer) C (/C£.; 2 (&t))* compactly for r > 0, and (/C£/ 2 (fcr))* C (£*/ 2 (fcT))* continuously 
for r < 1, by the Aubin-Simon Lemma (see e.g. [4], p. 98) with ([48]) and (|49|) . we obtain the 

relative compactness of the sequence ^°/ 2 (' CT )(/ »; Mi ;)) m L 2 {h, {^■ v a/2^ tT )T) f° r au r e ]0> ![• 
From (l2"5l) we deduce 



Mi* /2 (fcr)(A,«,) - P° /2 (kT)pu in L 2 (/ fe , (K£ /a (fcr))*) strong, Vr g]0, 1[. (50) 
Since from (|47|) , we have 

^°/ 2 (^)p;(t) = v* g h vr g]q, i[, (51) 

and we can write 

(p v u v ,Pv(t)(ur,))LZ(n)dt = / ( Pri u v ,P° /2 (kT)pr(t)u v } L 2 {a) dt 

( P a/2( kT )(Pri u v)' P *( t ) u v)L2(n)dt 

(P^ikrXp^^^u^^^^dt. 

The sequence (u v ) is bounded in L 2 (0,T, V) for all r G]0, 1[, therefore (P£(t)u n ) is bounded in 
L 2 (0,T,/C; /2 ) for all r G]0,1[. 

Therefore there exists a subsequence of P£(t)u v still denoted P£(t)u v such that 

i£(t)u, -» P£(t)u in i 2 (0,T,/C; /2 ) weak. (52) 
Passing to the limit in r/ yields 

^ Q (Pr,Ur,,P^(t)(Ur,)) L 2 in) dt = j {P° ^{kr) pu, F%(t)u) L '(0)dt 



( P U,P° /2 (kT)pr(t)u) L 2 {n) dt 

{pu,P T a {t)u) L 2 {n) dt. 



Summing over k = 1, .., JVr, we finally obtain 



which implies 



lim ||p,u, • P^K,) - • p;(m)|| l i (q) = 0, 

77— >U v ' 



lim lim • P£(u„) - p?i ■ P;(u)|| L i (Q) = 0. 

cr— »0 77—^0 



□ 



13 



We can now conclude to the strong convergence of u v . 
Let e > 0. From Lemma \TM 

3(T > 0, Vcr < (t , ||P;(u) - u\\ L 2 (Q) < e. 

From Lemma 13.71 

3a > 0, Vcr < ct , 3rj > 0, V77 < 770, ll-P^K,) - u»,||l 2 (q) < e, 
and by Lemma 13.81 

3t7 > 0, Vcr < cr , 3770 > 0, V77 < 770, Hp^u^P^u,,) - puP£(u)||z,i(g) < e. 
We therefore get from (up to the extraction of a subsequence) 

3i]o > 0, Vt] < 770, / \u n — u\ 2 dxdt < Ce 
Jo Jn 

which means that (still up to a subsequence) 

— > u in L 2 (Q) strong. (53) 

Classically, we also obtain from (fT9|) 

prjUrj — 1 p?j in L 2 (Q) weak. (54) 

3.4 Passing to the limit 

Let us now prove that as 77 goes to zero, a subsequence of (u Vl p v ) converges toward (u, p) solution 
of the weak formulation (fTT|) . Indeed : We have proved that p v — v p in i°°(0, T, L°°(r£)) weak * . 
Therefore 

P eL°°(Q). 

We have proved that u v — 1 tj in L 2 (0, T, 1/) weak, that ^p^u^ is bounded in L°°(0, T, L 2 (f2)) and 

p,, bounded from above and below in i°°(0, T, L°°(fi)). This implies that it,, bounded in L°°(0, T, £ 2 (fi)), 

thus its weak limit belongs to L°°(0, T, L 2 (Q)) n L 2 (0, T, V): 

I* e l°°(o, r, l 2 (^)) n i 2 (o, r, v). 

From Lemma 13.41 we have Hu = -ffit s = H(uq + lo x r) where H is the characteristic function of 
n s (t). Thus 

e K(t). 

Using compactness results of DiPerna-Lions we already obtained that p and H are solutions of 
transport equations with u and u s as velocities. For H this means that for all ip € C 1 (Q) with 
V>(T) = 0, 

/ / H^$- + Hu s ■ W efcrdt + / ffV(0) = 0. 
Jo 7n ot Jn 

As from Lemma [3~4l Hu s — Hu, H is also solution of 

/ \ H^- + Hu ■ Vip dxdt + [ H ip(Q) dx = 0. 
Jo Jn 9t J n 

In other terms H, like p satisfies a transport equation with velocity u. 
Let us finally check that u satisfies the momentum equation. 
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Let a > 0. If £t G H^Q) n L 2 (0, T; K\ (t)), from ([12]) and (JTSJ) we get 



dip u ) 1 

— + div(p,,u,, (g) u,,) - 2/idiv(D(u n )) + Vp ?) + -H v p v (u v - u v , s ) - p v g 
From Lemma 13.51 there exists 770 such that 77 < 770 implies: 

/ H v p ri (u v - u, hS ) ■£ (T dx = 0. 
Jo 

By integration by parts 

/ -2(j,div(D(u v )) -£ a dx = / 2pD(u v ) : D(£ a ) dx, 
Jn Jn 



■ £crdx = 0. 



/ div(p n u v (g) u n ) ■ £ CT dx = / -(p v u v ® u n ) : D(£ a ) dx, 
Jn Jn 



[ d(p v u n ) d f f d^ a 

In '^ dX= dti P ^ ^ adx ~ I P ^ ■ -dT ^ 



As a result 



Pr,u n ■ — + {p n u v ®u v ~ 2fiD(u n )) : D(£ a ) + p v g ■ 



dX = a\ / PvUri ' ^ a dX ' 



We have already established that u v u in L 2 (0, T, Hq (f2)) weak, u,, — > u in L 2 (0, T, L 2 (f2)) strong, 
p^w,, — 1 pu in L 2 (0, T, L 2 (Q)) weak, and p r) — > p in L 2 (0, T, i 2 (0)) strong. Letting 77 goes to zero, 
we thus obtain 

pu ■ — f + (pu(g)u- 2pD(u)) : D(£ a ) + pg ■ 
at 



dx = — pu - f CT da;. 



which corresponds to the weak formulation (fTTj) . This holds for any £ CT <E H 1 (Q) n £ 2 (0, T; if* (/;)), 
for arbitrary tr > 0, and, since the time interval has been chosen to guarantee that there is no 
contact with the boundary, by Proposition 4.3 of [19], for all £ G ^(Q) n L 2 (0,T; K(t)). This 
ends the proof of theorem 12.11 



4 Numerical simulations 

We give here a few numerical illustrations of the penalization method. We only sketch the nu- 
merical discretization and we refer to [3] for a more detailed description and further numerical 
results. 

We choose a time-step At and denote by a superscript n discretization of all quantities at time 
t n = nAt. For each time integration, we split the penalization model (fT2]) -(|16 p as follows 

• We solve the following variable density flow problem and obtain: 

u = u n - At(u n .V)u n + —pAu - —Vp n+1 + Atg 

p n p n 



with 



p n = Ps H n + p f (1-iP) 
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• We compute the rigid velocity the rigid velocity corresponding to u: 

p n uH n dx 



ug + w x r™ , uq 



p n H n dx 



UJ = J 



p n {r n x u)H n dx 



• We penalize the flow with this rigid velocity inside the solid body: 



At 



1 



H n (u s 



u n+1 ) & u n+1 



n 



i + 



• We finally advect the solid with the rigid velocity: 

H n+1 = H n -Atu„-VH r ' 



Note that we have used an implicit time discretization of the velocity penalization, which allows 
to use a very small penalization parameter rj. Using an explicit method would require this value 
no to be smaller than At. It can indeed be checked that a explicit method with rj = At essentially 
amounts to the projection method [17] • We will see below that using smaller values of rj together 
with an implicit scheme has a significant effect on the accuracy of the method. 

In order to numerically validate the penalization method, we consider the case of the sedimen- 
tation of a rigid cylinder in two dimensions (see pj], [7]). The domain fl — [0,2] x [0,6] is filled 
with an incompressible viscous fluid initially at rest, of density pf = 1 and viscosity jXf — 0.01. 
The rigid cylinder of radius R — 0.125 and density p s — 1.5 is initially centered in (1,4), and we 
apply the gravity force g — —980. 

In order to verify how the rigid constraint is satisfied in the solid, we monitor at time t = 0.1 the 
L 2 -norm of the discrete deformation tensor defined by: 



\\D{u)f 



ij,4>ij<a 



-12 



at 



We fix Ax — 1/256 and At = 10 4 , and compute this norm for values of 77 from 10 4 to 10 
t = 0.1. 

The results presented in table [T] indicate a convergence of the method with order 1 in r\. 
In figure [1] we show the profiles of the vertical velocity for several values of 77, corresponding to a 



V 


\\D(u)\\ L 2 {ns{t)) 


a for 0(-q a ) 


10- 


-4 


4.30838 




10" 


-6 


3.84749 x 10- 2 


1.0247 


10" 


-8 


3.45379 x 10~ 4 


1.0234 


10- 


10 


3.81643 x 10~ 6 


0.9783 


10- 


12 


3.79832 x 10~ 8 


1.001 



Table 1: Sedimentation of a two dimensional cylinder. Errors on ||D(u)||£2(Q s ( t )) and convergence 
orders at t = 0.1 for Ax = 1/256 and At = lO" 4 . 



cross section at the center of the cylinder. We can observe that below rj = 10 8 one may consider 
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that we obtained converged velocity results. By taking 77 = 10 -4 , that is rj — At, along with an 
explicit treatment of the penalization term, we get the projection method. As far as precision 
is concerned, one can note the benefit of using larger penalization parameters combined with an 
implicit time discretization of the penalization term. 




-12 - 



0.5 1 1.5 2 

x 



Figure 1: Sedimentation of a two-dimensional cylinder, for Ax = 1/256 and At = 1CP 4 . Vertical 
velocity in an horizontal cross-section through the center of the cylinder at t = 0.1 for several 
values of the penalization parameter. 



5 Conclusion 

We have presented and analyzed a penalization method that extends the method of [5] to the 
case of a rigid body moving freely in an incompressible fluid. The proof is based on compactness 
arguments. Numerical illustrations have been provided to illustrate our convergence result. The 
benefit of using very large penalization parameters combined with an implicit time discretization 
of the penalization term, compared to the projection method |21j which corresponds to a particular 
explicit time discretization for the penalized equation, has been demonstrated. 

While this was not our primary goal, an outcome of our convergence study is an existence 
result for a weak formulation of the coupling between a rigid solid and a fluid. Let us shortly 
discuss how this result compares with existing ones [12l [8] [6] [19] . In [12] local in time existence 
and uniqueness of strong solutions was proved. The Eulerian approach was developed in [5] where 
global in time existence of weak solutions was proved in dimension 2, without collisions. In 
the three-dimensional case, to our knowledge only local in time existence of weak solutions was 
obtained, since L 2 regularity of the time derivative of velocity was required (and therefore global 
existence would imply global existence of strong solutions) . In [6] [13] the existence of global weak 
solution in three dimensions for one ball shaped solid, with possible collision with the boundary, 
was proved. In [19] the existence of global weak solutions for several rigid bodies with collisions 
was proved in dimension 2. Our results prove the existence of global in time weak solutions in three 
dimensions, before collision. By contrast with J5][T3], this result can easily be generalized to the 
case of several bodies by introducing indicator functions, rigid velocities and penalization terms 
corresponding to each body. To our knowledge, the existence of global in time weak solutions for 
several bodies with collisions is an open problem in three dimensions. 

6 Appendix 

This section si devoted to the proof of some technical lemmas that were used in section 3. 
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Proof of Lemma 3.3 Let <p(t) S L 2 (f2\fi Sl0 -(t)), and (v(t), q(t)) e H 1 (Q\n a<a (t))xL 2 (Q\Q Sl(T (t)) 
solution of the Stokes problem 



-Av(t) + Vq(t) = <f>(t) on Q\VL s ^{t) 
divv(f) = 

K*) = ° 



on f2\n SiCT (t), 
on 



9 (n\n, )0 .(t)) . 

Since we assumed a C 2 regularity on f2\fig this regularity is conserved through rigid motion and, 
for a small enough (say a < <7o for some Co > 0) to n\fl s . a (t). The regularity results of Agmon- 
Douglis-Nirenberg on the linear Stokes problem (see 22J, prop. 2.3. p. 35) give 

(«(t),g(t)) g ff 2 (r>\n77(t)) x fl- x (n\rv(t)) 

and there exists C > such that 

+ ll?(*)lli?i(Q\!v7(t)) ^ c ll#)llz,=(o.\T^(i))' 

Note that, with our definition of T, the constant C, which depends on the geometry of the domain 
boundary, can be taken independent of t G [0, T] and er, provided do is taken small enough. We 
can then write 



/ w{t) ■ cj)(t) dx = - / w(t) • Au(t) da; + / w(t)-X7q(t) 

Ja.\n s ,v(t) Jn\n.,_„(t) Jn\n.,.„(t) 



dx 



i(l\(i,, a (t) 



a(f2\n s ,„)(t) 



n\n i><r (t) 



, s dv(t) , 



n\n, 1<r (t) 



Vw(i) • Vu(t) da; 



9(o\o s ,„)(t) 



w(t)q(t) ■ nds 



a\n,, a (t) 



div w(t)q(t) dx 



, , 9u(t) , 
t«(t) • -^ds 
an 



9(o\a s ,„)(t) 



dn 



v(t)ds 



n\n,,<r(i) 



9f2„,<,(t) 



w(t)q(t) ■ nds 



an a , CT (t) 



an 3 ,„(t) 



w(t)q(t) ■ nds 



v(t) ■ Vp(t) dx 



The integral of i> • Vp vanishes since i> is divergence-free and vanishes on d (Q\n SiCT (t)). Then 
using classical trace theorems in Sobolev spaces we get 



/_ 



w(i) • </>(t) da; 



an s ,„(t) 



5 (t)--^d s 



g(t)q(t) ■ nds 



an a ,<,(t) 



< HsWIU^an.,^*)) (||Vw(i)|| L 2 (a0s ct(4)) + ||g(*)||i2(an iltr (t))) 

< C||g(i)|| L 2 (ans ^t)) (jK£)|| ff 2 (fA f^ (t)) + lk(*)|| H i(n\fi^:(t)) 

< C||5(*)IU 2 (sn,,„(t)) II <H*) llL2 (0 \n77(t)) 

1 1 

< c 'll«(*)III»(n. 1<r {t))ll Vtt (*)IIE«(n., ir (t))ll^( t )ll£»(n\I5T(t))- 



This proves the assertion. 
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Proof of Lemma 3.6 

Step 1: We first show how to construct for a.e. t G]0, T[ a function u CT (.,i) £ /C£(i) such that 

lim - u(.,t)||v = a.e. on ]0,T[. 

Let cr > and w CT (., i) such that 



(-Av a (.,t) + Vp(.,t) 
div t) = 
u CT (.,t) = u s (.,t) 

K(.,t) = o 



-Au(.,t) on n\n Sj(7 (t), 
on n\n SjCT (i), 

on <90 S)O .(t), 
on <9J7, 



where 



p{r x da; I x r. 



By lemma I3T41 i) = u s (., t) on fi s (t). Extending u CT (., t) by u s (., t) in ^s,a(t), we have i><r(., £) € 
JC r a {t). We set eo-(.,i) = i;<r(.,t) — u(.,t). It satisfies 



'-Ae a (.,t) + Vp(.,i) = on n\n,, a (t), 

dive CT (.,t) = on f2\f2 S)(T (£), 

e CT (., i) = u s (., i) - i) on 9f2 8)<7 (*), 

,e (7 (. ! <) = on 90. 

We extend e CT (.,i) by u s (.,t) — u(.,t) in fi Si<T (t), so that e a (.,t) = in fi s (i). 
We claim that 

lim ||e (T (.,t)|| L 2 (!:2) = a.e. on ]0,T[. 



(55) 



In tt s (t) e a (.,t) = 0, thus 

ll e <r(-^)|| i 2( ) = H e cr(-^)llL 2 (f2 s ,„(t)\n7(t)) + H e o- *)H i 2 (f7\Ol7T(*)) - 

Since Q Sj(7 (t)\Q s (t) has width 2a, from the proof of lemma 5.10 of [10] we have a.e. on ]0, T[, 



\K(;t)\\ 



"i/ ii w (ll e c r (-,t)|| L 2( 9f2s ( t ) ) + o- 2 ||Ve CT (., i)|| L 2( 0s (T(t) \ ns(t )) 



Ci 



c 



ll e <T(-,*)|| L 2( S:23 ( t ))||Ve (J (., i)|| L 2( ns ( t )) + cr 2 ||Ve a (., t)|| i2( ^ s ,„(t)\n s (t))) 



Ca 2 ||Ve CT (.,i)|| 



L2(f2 s ,„(t)\f2 s (t)) 



Next, as e a (.. t) = u s (., t)—u(., t) in f2 S]Cr (t) and t) and it s (., i) are in Hq (f2), we get ||Ve CT ( 
||Ve CT (., f)|| i 2( f2s / f \\ < || Ve CT (., i)||j,2(m < C, where C is independent of cr. This gives 

ll e <x(-, *)lll,2(n S)<T (t)) = ll e cr(., *) lli,2(0 Sj(r (t)\m(t)) - ( ^ Cr ' 

By Lemma l3~3l we thus get 

ll e ^(-^)lli2 (fi \f^ (t)) <C'||e (T (.,i)|| i 2 (0si(TW) ||Ve (J (.,<)|| L 2 (f2s ff(t)) < Ccr. 
Collecting the above estimates, we conclude that 

lim \\e a (.,t)f L2in) = 0. 

In order to prove that this convergence also holds in V r we first note that 

||e CT (.,t)|| H1(fi) <Ca.e. on]0,T[ (56) 



li 2 (n., ff (t)\n.(t)) 



< 
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as is readily seen from estimates on the Stokes problem verified by e a . By interpolation (see e.g. 
PP, p. f35), we obtain 

||e„(.,t)||vr < ||e CT (.,t)|| (57) 

and due to ([55*)) and (f5"""[) , 

lim ||e CT (.,f)|| V r = Vr G [Q,l[ a.e. on]0,T[. (58) 
Step 2: By definition of P£, 

\\pr U (.,t)-u(.,t)\\ Vr < \\v a (;t) - U(.,t)\\ vr 

thus the pointwise convergence on v a we just obtained implies 

lim ||P>(.,f) - u(.,t)\\ V r = a.e. on ]0,T[. 

Step 3: f a : 1 1— > ||P£n(.,i) — u(.,t)||yr is measurable on [0,T] and since G /C£, 

l|p; u (,t)-<,t)ll| r di < / h(.,t)|itdt 



(59) 



II/ct|I lI(o,t) : J 



2(l~r-) 



< c IK.,t)|| £B f n) |K,i)ll 



dt 



< C|| M llLcxr( 0jTiL 2( n ))||u|| L 2( 0!T!fl -l( n )) 

c. 



To summarize / CT verifies 



lim Ut) = a.e. on [0,T], 

f a is measurable on [0, T], 
ll/JI 2 < C with r < 1. 

l|J5|l If (0,T) — 

Therefore, thanks to lemma [3~"2l lim ll/o-llraco t) = 0' which means 

lim \\Py - u\\ L 2 {0 ^ T ^ r) = 0. 



Proof of Lemma 3.7 

Step 1: We construct for a.e. fixed t € [0, T] a function v v<J (.,t) £ K, r a (t) such that 
lim lim ||iw(., t) — u. n {., t)||y = a.e. on ]0, T[. 

cr— >0 r\— >0 

Let g > and v V(7 (.,t) solution of the following Stokes problem outside f2^ CT (i): 



-Atv(.,t)+Vp(.,i) = -Au,(.,t) on 
divu, 7CT (.,t) = on Q\SfH ;(T (t), 

v v<r (.,t) =u ntS (.,t) ondQ.^ (t), 



^v va (.,t) = 



on dft. 



Extending by u VjS (.,t) in ^^(i), we have v v<T (.,t) € IC a (t). We then introduce e va (.,t) = 

v r]a (.,t) — u v (.,t). It verifies 



-Ae,^(.,i)+Vp(.,t)=0 (mfl\n^(i), 
dive J)(T (.,t) = on fi\fil? i0 .(t), 

e I)cr (.,i) = u, hS (.,t) ~ u v (.,t) on df^^i), 



e I)(J (.,i) = 



0fi, 
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and we extend it by u^ )S (.,t) — u v (.,t) in ffli a (t). 
We claim that 

lim lim \\er,*(.,t)\\i*(n) = a.e. on]0,T[. (60) 
From lemma 15751 for a given a > 0, there exists 7/0 > such that V77 < 770, 

nj(t) c n. )ff (t) and n»(t) c fi2 i<y (t). 

Let r\ < i]o. We write 

|| e T7C r(-, i)||i2(n) = ll ej / CT (-' *)lli 2 (nj(t)) + ll e w(-' i )llL2(n s-(T (t)\o7(t)) + II^C *)II.L 2 (n\n^7(t))- (61) 
From estimate there holds 

Il e ^(-»*)llia(n2(t)) d * ^ ^ ( 62 ) 
Since r2 s CT (t)\fi2(i) has width less than 2cr, from the proof of lemma 5.10 of [TD] we have a.e. on 

]o,n 

ll e »?<T(-^)llL2(o 3 ^(t)\?27(t)) — ^ (ll e ^(-' *)Hx,2(af22(t)) + fj2 ll^ 7e ')CT(-: *)Hi2(n SiCr (t)\nJ(t))) ■ 

(63) 

And using a trace theorem, we get 

1 1 e^cr (. , £)|| L 2(q s ier (i)\nf(t)) - ^ (ll e w(-' i )llL2( j( t ))||Ve, )cr (.,t)|| i2(fi 2(t)) 

+o- 2 ||Ve ?)cr (., i)|| L 2(n s l0 .(t)\nT(t))) 

< C (ll e »)cr(-^)|| L 2( ' s ;( t )) ||Ve ?)cr (., *)llx^(n, i<r (t)) 

+a 2 ||Ve w (.,i)|| 2 L2(0s ^ (()) ). (64) 

Adding ||e I)cr (., i)|| 2 to this inequality gives 

ll e »)cr(-^)llL 2 (n s ,„(t)) - ^ (ll er ; CT (-' *)Hz, 2 (n3(t)) + ll e ^(-' i )llL 2 (oj(t))ll^ 7ei ; cr (-' <: )llL 2 (n Si „(t)) 

+a 2 ||V e , CT (.,t)|| 2 2(ns ^ (t)) ). (65) 



For the last term in (|6ip we use Lemma 

ll e »)o-(-^)llL 2 (J2\f^7(i)) - ( ^ll e '? CT (-' i )llL 2 (n s , <T (t))ll^ 7e 'J cr (-'*)llL 2 (fi 3 .„(t))- (^) 

Since e. (?CT (., i) = u^ )S (., i) — u^(-, i) in fi S)(T (i) and u^, u VtS are bounded in L 2 (0, T; iJp (f2)), we have 

T ||V e , CT (.,i)||i 2(0 _ (t)) *<C. (67) 

With (|6"2j) and (|6"T)) we are now in position to estimate the integral over [0,T] of (|64j|66|) . By 
Cauchy-Schwarz inequality: 

T 
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J o \M,mh { n\n^(t)) dt ^ C yJ o \M,t)\\Un s , At)) dtj <C( V + ^+a^. 

Therefore, for a fixed value of a we can pass to the limit in r\, and then pass to the limit in a , to 
obtain 

rT 

IM->*)l&(n)* = °- ( 68 ) 



lim lim 

cr->0 r;-»0 



This strong convergence can be turned into an almost everywhere in t convergence up to the 
extraction of a subsequence. The rest of the proof is adapted in a straightforward way from that 
of Lemma [ 



Acknowledgments 

This work was supported by the French Ministry of Education through ANR grant 06-BLAN-0306. 

References 

[1] R.A. Adams and J. Fournier, Sobolev spaces, second edition, Elsevier (2003) 

[2] Ph. Angot, C.-H. Bruneau and P. Fabrie, A penalization method to take into account obstacles 
in incompressible viscous flows, Numer. Math. 81: 497-520 (1999) 

[3] C. Bost, Methodes Level-Set et penalisation pour le calcul d 'interactions fluide- structure, PhD 
Thesis, University of Grenoble, France (2008) 

[4] F. Boyer and P. Fabrie, Elements d 'analyse pour V etude de quelques modeles d'ecoulements 
de fluides visqueux incompressibles, Mathematiques & Applications 52 (2006) 

[5] H. Brezis. Analyse fonctionnelle, Theorie et applications, Masson (1992). 

[6] C. Conca, H.J. San Martin and M. Tucsnak, Existence of solutions for the equations modelling 
the motion of a rigid body in a viscous fluid, Comm. Partial Differential Equations 25, 1019- 
1042 (2000) 

[7] M. Coquerelle and G.-H. Cottet, A vortex level-set method for the two-way coupling of an 
incompressible fluid with colliding rigid bodies, Journal of Computational Physics, 227, 9121- 
9137 (2008) 

[8] B. Desjardins and M.J. Esteban, Existence of Weak Solutions for the Motion of Rigid Bodies 
in a Viscous Fluid, Arch. Rational Mech. Anal. 146, 59-71 (1999) 

[9] R.J. Di Perna and P-L. Lions, Ordinary differential equations, transport theory and Sobolev 
spaces, Invent, math. 98, 511-547 (1989) 

[10] H. Fujita and N. Sauer, On existence of weak solutions of the Navier-Stokes equations in 
regions with moving boundaries, J. Fac. Sci. Univ. Tokyo Sec. 1A (became from 1993 Journal 
of mathematical sciences, the University of Tokyo), 17, 403-420 (1970) 

[11] R. Glowinski, T. W. Pan, T. I. Hesla, D. D. Joseph, and J. Periaux, A Fictitious Domain 
Approach to the Direct Numerical Simulation of Incompressible Viscous Flow past Moving 
Rigid Bodies: Application to Particulate Flow, Journal of Computational Physics 169, 363- 
426 (2001) 

[12] C. Grandmont and Y. Maday, Existence de solutions d'un probleme de couplage fluide- 
structure bidimensionel instationnaire, C. R. Acad. Sci Paris Ser. I Math. 326, 525-530 (1998) 



22 



[13] M.D. Gunzburgcr, H.-C. Lee and G.A. Seregin, Global Existence of Weak Solutions for Viscous 
Incompressible Flows around a Moving Rigid Body in Three Dimensions, J. math, fluid mech. 
2, 219-266 (2000) 

[14] M. Hillairet, Lack of collision between solid bodies in a 2D constant density incompressible 
viscous flow, Communications in Partial Differential Equations, 32:9, 1345-1371 (2007) 

[15] J. Janela, A. Lefebvre and B. Maury, A Penalty Method for the simulation of fluid-body rigid 
interaction, ESAIM: Proceedings, Vol.14, 115-123 (2005) 

[16] J.-L. Lions, Quelques methodes de resolutions des problemes aux limites non Uneaires, Dunod 
(1968) 

[17] N.A. Patankar, A formulation for fast computations of rigid particulate flows, Center Turbul. 
Res., Ann. Res. Briefs, 185-196 (2001) 

[18] J. San Martm, J.-F. Schcid, T. Takahashi, and M. Tucsnak, Convergence of the Lagrange- 
Galerkin method for the equations modelling the motion of a fluid-rigid system, SIAM J. 
Numcr. Anal., 43 (2005), pp. 1536-1571. 

[19] J. A. San Martin, V. Starovoitov and M. Tucsnak. Global Weak Solutions for the Two- 
Dimensional Motion of Several Rigid Bodies in an Incompressible Viscous Fluid, Arch. Ra- 
tional Mech. Anal. 161, 113-147 (2002) 

[20] A. Sarthou, S. Vincent, JP. Caltagirone and Ph. Angot, Eulerian-Lagrangian grid coupling 
and penalty methods for the simulation of multiphase flows interacting with complex objects, 
Int. J. Numer. Meth. Fluids 00:1-6 (2007) 

[21] N. Sharma and N.A. Patankar, A fast computation technique for the direct numerical simu- 
lation of rigid particulate flows, Journal of Computational Physics 205, 439-457 (2005) 

[22] R. Temam, Navier-Stokes equations and numerical analysis, North-Holland, Amsterdam 
(1979) 



23 



